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SUMMARY 

The NASA Lewis heat pipe code has been developed to predict the performance of heat pipes in 
the steady state. The code can be used as a design tool on a personal computer or, with a suitable 
calling routine, as a subroutine for a mainframe radiator code. A variety of wick structures, including 
a user input option, can be used. Heat pipes with multiple evaporators, condensers, and adiabatic sec- 
tions in series and with wick structures that differ among sections can be modeled. Several working 
fluids can be chosen, including potassium, sodium, and lithium, for which the monomer-dimer equili- 
brium is considered. The code incorporates a vapor flow algorithm that treats compressibility and 
axially varying heat input. This code facilitates the determination of heat pipe operating temperatures 
and heat pipe limits that may be encountered at the specified heat input and environment temperature. 
Data are input to the computer through a user-interactive input subroutine. Output, such as liquid and 
vapor pressures and temperatures, is printed at equally spaced axial positions along the pipe as deter- 
mined by the user. 


INTRODUCTION 

The NASA Lewis Research Center is conducting research programs involving power production 
for space applications. Systems being considered may employ heat pipes either to conduct heat from 
the source to the power converter or to remove waste heat to radiators. It was desired to extend the 
capabilities of current steady-state codes (refs. 1 to 3) in order to perform the design and analysis 
required for evaluating systems of interest. Hence new computation algorithms were derived and a 
new code was written that incorporates portions of an existing Lewis code (ref. 3). This code has 
features that should enhance its execution speed and precision. 

An important requirement was a vapor flow algorithm that could treat compressibility and axially 
varying heat input in a satisfactory manner yet be concise enough to function on a microcomputer of 
limited capacity and speed. Finite difference techniques, which had been used effectively in funda- 
mental studies of heat pipe vapor dynamics (refs. 4 and 5), were eliminated by the latter consideration. 
Also desired was great flexibility in the input of boundary conditions and other parameters, so that a 
wide variety of cases could be computed. We believe that for our purposes this code represents an 
improvement over other similar codes. 



A listing of the NASA Lewis heat pipe code is available from the Computer Software Manage- 
ment and Information Center, Athens, Georgia 30602. Because of possible updates of the code, the 
results of example cases run using the listing code may differ slightly from those presented herein. 

The authors are indebted to Barbara McKissock and Carlos Rodriguez for their help and patience 
in running the code during the debugging process. 


DISCUSSION 

Guidelines established for the Lewis code were, among others, a simple, logical, modular structure 
and extensive documentation. Also, a code was sought that could be used as a design tool on a per- 
sonal computer (PC) or, with a suitable calling routine, as a subroutine for a mainframe radiator code. 

The Lewis code currently has about 45 subroutines, each providing a limited function, that are 
controlled by a few master calling routines. This structure makes program flow easy to interpret. 
Moreover, it facilitates the replacement or updating of subroutines and the addition of new ones. Most 
of the required physical parameters, wick resistance for example, are computed from input specifica- 
tions for particular types of structures by using accepted formulas, as in other codes. For the benefit 
of the user who wishes to consider internal structures outside of the repertory included, provision is 
made for user input of the required physical parameters. 

The general features of the code are described first, followed by a brief discussion of its func- 
tioning. Appendix A presents a list of symbols appearing in the input and output of the code. 
Appendix B contains a discussion of the vapor flow algorithm employed. Appendix C describes 
and illustrates the user-interactive data entry for an example case, while appendix D describes the 
data output for the example. Appendix E gives the output for another example case. Appendix F 
summarizes the significance of the various heat pipe limits. The code was written in Fortran 77 and 
compiled with Microsoft Fortran 4.1. 


Summary of Code Features 

The code has been constructed for straight, cylindrical heat pipes with no bends. Simple modi- 
fication to noncircular cross sections could be achieved by using the hydraulic radius approach, but 
this feature is not currently implemented. Curvature, but no sharp bends, is permissible. Operation in 
the presence of gravitational fields is accounted for in the case of straight pipes, but thermosyphons 
are not included. 

The heat pipe can be subdivided into as many as 20 sections, far more than will normally be 
required. These sections can be heat input or evaporator sections, adiabatic sections, and heat removal 
or condenser sections, mixed in any manner provided that the first section is a heat input section. The 
nature of any section is indicated by the type of boundary condition data that is entered for that 
section. 

At present the following options are available for the boundary conditions: (1) specified heat input 
or removal rate at the pipe exterior surface and (2) specified environment temperature external to the 
pipe, with heat transfer conditions from the pipe surface to the environment provided by the user. 
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For each section the required data can be entered as single values pertaining to the entire section. 
Values of the input parameters can also be entered at as many as 20 points spaced along the section, 
with a different magnitude at each point. 

Computations along the pipe are made by a Runge-Kutta routine for which the initial step size is 
chosen, with automatic reduction or increase as the slope of the principal dependent variable, pressure, 
exceeds or falls below certain bounds. 

Data input is by means of an interactive subroutine that queries the user concerning the options to 
be employed in the case to be run. For output the user can select the number of uniformly spaced 
points along the pipe axis at which the principal variables will be printed out. The spacing of these 
points need not correspond with the step size chosen for the Runge-Kutta solution. The printout is 
furnished after thermal convergence of the solution has been obtained. When a heat pipe limit is 
encountered during the solution that does not cause the solution to halt, the type of limit is printed out 
with the solution. The data printed out can be output directly to a plotting device if a suitable plotting 
utility is available and the necessary changes are made in the source code. 

Subdivisions of Heat Pipe 

To clarify the various types of subdivisions of the heat pipe, consider this illustration of a heat 
pipe, with two evaporators separated by an adiabatic section, and a condenser: 

Pipe sections . — The major subdivisions are the pipe sections, which can be evaporator, adiabatic 
section, and condenser in the usual terminology. So that other than the usual sequence of boundary 
conditions can be employed, this nomenclature is not used in the code. Rather, the pipe is conceived 
to be divided into NS sections numbered IS = 1 to NS. The function of each section is denoted by the 
type of boundary condition imposed at the time of data input. The number of sections is limited to 
20 . 


Uniform boundary conditions pertaining to the entire section can be entered at the time of data 
input. Alternatively, varying conditions can be entered at I locations along each section numbered 
IS. The total of such locations for each section is designated by NH(IS), where NH(IS) must be less 
than 21. At any section end that abuts another section, the same physical coordinate can have two 
boundary conditions associated with it: one pertaining to the section ending, the other to the next 
section. 

Runge-Kutta steps . — The initial number of steps, ISTEP, that are employed in the Runge-Kutta 
solution is chosen by the user. This choice determines the maximum step size to be employed in the 
program. When the pressure drop between two successive steps exceeds about 3 percent, the subse- 
quent step size is subject to automatic halving during execution. This is particularly important when 
conditions change rapidly as choking is approached. Similarly if the pressure increase is less than 1 
percent, the step size is doubled unless the maximum step size has been reached. Also, when the end 
of a section is approached, the step size is adjusted so that a step coincides with the end of each 
section. 

Print locations . — For purposes of data output the user can choose the number of equally spaced 
locations, numbered IP, along the pipe at which data will be printed out after a solution has been 
obtained. The total number of such locations is IPSTEP + 1, where IPSTEP is less than 101. For all 
locations at which a print location and a Runga-Kutta step do not coincide, the data for printing are 
linearly interpolated between the adjacent Runge-Kutta steps. Also printed out in the proper axial 
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sequence are the data pertaining to the beginning and end of each section whether or not these 
correspond with a chosen print location. Except for the pipe beginning and end, there will be two 
listings of data with the same axial coordinate at each such section change. 


Boundary Conditions in Sections 

Each section of the pipe can have imposed on it one of three boundary conditions: 

(1) specified heat input at the exterior surface, (2) temperature at the exterior surface, and (3) 
environment source or sink temperature. Options 1 and 3 are fully implemented at this time, but 
option 2 is not. 

Option 1, specified heat input, requires information on heat input. For sections with uniform radial 
heat input, only the total heat in watts for the section is required. If the radial heat input varies along 
the section, the input rate in watts per meter at as many as 20 points along the section can be entered. 
The cubic spline interpolation produces a smooth fit of the data. A section with no heat input will be 
an adiabatic section. 

Option 3, specified environment source or sink temperature, requires that parameters pertaining to 
heat transfer between the pipe surface and the environment be provided. In the case of radiation these 
are the specific area of the pipe surface in square meters per meter of length and a radiation constant 
that combines the view factor and the effective emissivity. In the case of convection or conduction a 
suitable heat transfer coefficient must be provided. Combined convection and radiation in a section 
can be handled. Again these data can be entered as single values for an entire section. If there is 
variation along the section the data can be entered at several points, with spline interpolation being 
used to provide a smooth fit between the points. 

A pipe can have sections with both options for heat input or both options for heat removal. For 
example, a section with prescribed heat input and a section with heat gain from a source by radiation 
can be present in the same pipe. 

The user should be aware that not all conditions can be handled. There is an obvious requirement 
that any heat pipe have at least one heat input section and one heat removal section. Furthermore 
some choices of boundary conditions may result in an encounter with one of the heat pipe limits, 
necessitating a change in the dimensions or the working fluid. 


Spline Interpolation of Input Data 

The flexibility of the program is enhanced greatly by the use of cubic spline interpolation of some 
of the input data. With this approximation, data that are expected to be continuous but varying with 
position in a section need be entered for only a few points. Coefficients are then determined that 
enable the data to be interpolated smoothly at any axial location within the section. Among the data 
that use the spline fit are the radiation parameters and convective heat transfer coefficients for the pipe 
surface. Specified heat input that varies with position in the section can also be interpolated by the 
spline fit provided that the change within the section is smooth. The algorithms to implement the 
spline interpolation were adapted from reference 6. 
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Material and Working Fluid Choices 


The code contains thermal conductivity data for the following wall and wick materials: aluminum, 
beryllium, copper, Inconel, iron, molybdenum, nickel, niobium-1 % zirconium, 304 stainless steel, tan- 
talum, titanium, and tungsten. If the material of interest is not in this list, the thermal conductivities 
of the wall and wick materials can be entered by the user from other sources. Similarly the effective 
thermal conductivity of the fluid-saturated wick can be entered in the user-input wick property section 
of the code. The properties so entered will be for some estimated mean operating temperature of the 
pipe and will be subject to correction. 

The following working fluids are provided in the code: alkali metals (lithium, potassium, sodium, 
and cesium) and other fluids (ammonia, methanol, mercury, Dowtherm A, acetone, benzene, Freon- 1 1, 
n-heptane, silver, toluene, and m-xylene). 

Thermodynamic Data for Working Fluids 

Alkali liquid metals. — The alkali liquid-metal vapor contains monatomic (monomer), diatomic 
(dimer), and tetrameric species (ref. 7) . The Lewis heat pipe code requires that the monomer and 
dimer species be accounted for in computing the composition and thermal properties of the working 
fluid. 

The required thermodynamic data for the alkali metals were taken from reference 8. Reference 8 
presents polynomial least-squares Fits for the specific heats of the solid, liquid, and gaseous species as 
a function of temperature. The data for the gaseous species pertain to the ideal-gas state of the indivi- 
dual species. Reference 8 also lists the enthalpies and entropies of the species relative to the atomic 
constituents, assigned at a standard temperature. The polynomial coefficients and the reference enthal- 
pies and entropies are stored in the heat pipe code as block data 

By using the polynomial coefficients and the assigned thermal data, the specific heat, enthalpy, and 
Gibbs free energy of each species can be computed as a function of temperature. A call to one of the 
thermal data subroutines at a specified temperature determines such properties as the composition of 
the vapor, the heat of vaporization, the vapor pressure, the speed of sound, and several first partial 
derivatives used in the code. 

In computing the properties the assumption is made that the vapor behaves as a mixture of ideal, 
but reacting, gases. This approach is similar to one described in reference 9 and referred to as the 
“quasi-chemical equation of state.” The assumption was made in reference 9 that monomeric, dimeric, 
and tetrameric species were present as ideal gases, but the data were fitted to an experimentally 
determined PVT diagram. In the present code the source of the thermal data is a tabulation of proper- 
ties for gaseous species that is based upon spectroscopic data. The tetramers, present in only a small 
fraction of the dimer concentration, are neglected. Computations of vapor pressure and heat of vapor- 
ization that are made by using the algorithm of the Lewis code agree with tabulations from other 
sources (refs. 10 and 1 1). Agreement is within 4 percent, which is within experimental error for data 
taken on heat pipes. 

Thermodynamic dat a for other working fluids . — Thermodynamic data for the other fluids are from 
references 10 and 11. 
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Physical Property Data 

The conductivities of the solid materials and the physical properties of the working fluids are from 
references 10 to 15. 


Flow Algorithms 

Vapor flow algorithm . — The vapor flow algorithm employed in the code was rederived for com- 
pressible flow from expressions for incompressible flow by Busse (ref. 16). Appendix B contains the 
derivation presented in NASA CR-185179 (ref. 17), with additions and corrections. 

For alkali metal working fluids, chemical equilibrium exists among the monomeric and dimeric 
gaseous species. No tangential velocity component exists at the liquid-vapor interface, and the radial 
pressure gradient may be neglected. 

The incompressible Busse equations (ref. 16) have been rederived as compressible relations for use 
in the present code as shown in appendix B. The change in pressure and the Busse velocity profile 
factor, which adjusts the Poiseuille profile, are computed at each step with density variation consid- 
ered. In addition, for alkali liquid metals the mean gas temperature, which is different from the 
saturated vapor-liquid temperature, is also determined at each step, by means of an energy integral 
approach. For these working fluids, large temperature differences can exist between the liquid-vapor 
interface and the core of the vapor (ref. 5). However, for these fluids the mean vapor temperature and 
the vapor-liquid interface temperature are equated beyond the point in the condenser where they 
become equal. For all working fluids at peak Mach numbers less than 0.01, Busse’s incompressible 
relations are employed, and the density at the end of each Runge-Kutta step is updated in a procedure 
similar to that employed in reference 1 for all Mach numbers. In this way a mathematical problem 
was avoided at conditions where the change in pressure and density is extremely small. 

At this time the assumption is made of laminar flow though the pipe. Studies (ref. 1 8) have 
shown that in regions of vapor injection, as in the evaporator, turbulence is suppressed at axial Rey- 
nolds numbers of up to 1 million. In regions of vapor extraction (condenser), turbulence is induced. 

A reliable theory as to the axial rate at which turbulence develops would be required for this effect to 
be incorporated into the code. Many other uncertainties concerning the condenser exist, such as recir- 
culation and droplet nucleation. An attempt to account for turbulence would seem premature until suf- 
ficient experimental data are available to clarify this region better. 

Liquid flow algorithm . — Liquid pressure drop is computed by means of the familiar Darcy equa- 
tion (ref. 10). Multiple axial liquid paths of differing permeability can be accounted for, as in artery 
structures. Hydrostatic pressure changes due to tilting of the pipe in gravity fields are considered. 


Wick Structure Options 

There are six wick options in the code: user input, screen wick, screen wick with arteries, sintered 
metal, sintered metal with arteries, and longitudinal grooves (not implemented). 

The user must select one of these options. Figure 2 shows the wick structures that are incorpo- 
rated in the code. The user input option, which is not shown in figure 2, has been incorporated to 
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allow the user to analyze heat pipes with wick structures other than those listed. The user input option 
requires the following data: 

(1) The wick permeability to fluid flow multiplied by the wick cross-sectional area 
perpendicular to fluid flow 

(2) The perimeter of the wick at the liquid-vapor interface 

(3) The porosity of the wick 

(4) The effective thermal conductivity of the wick 

(5) The effective vapor area perpendicular to the axis 

(6) The minimum effective capillary radius 

(7) The entrainment critical dimension 

(8) The minimum nucleation radius 

Option 2 is for screen wick structures as shown in figure 2. The screen wick option requires the 
following input data: 

(1) The diameter of the screen wire 

(2) The width of space between screen wires 

(3) The crimping factor 

(4) The number of screen wraps 

Option 3 can be used for screen wicks with arteries. This case requires the data from option 2 
plus the following: 

(1) The artery external diameter 

(2) The artery wall thickness 

(3) The number of arteries 

Option 4 can be used to analyze heat pipes with sintered-metal wick structures. This option 
requires the following data: 

(1) The porosity of the wick (the volume of liquid in the wick divided by the total volume of the 
wick) 

(2) The effective radius of the wick particles 

(3) The thermal conductivity experimental constant, which is 0.34 for metal felt or fiber and 0.53 
for sintered powders. The experiments and the derivation of the thermal conductivity equation 
are outlined in reference 10. 

(4) The thickness of the wick 

Option 5 can be used for sintered wicks with arteries. This case requires the data from option 4 
plus the following: 

(1) The artery inside diameter 

(2) The number of arteries 

(3) The protrusion area of each artery into the vapor space, which is used to account 
for the reduction in vapor flow area caused by the arteries as shown in figure 2 

Option 6 is for heat pipes with longitudinal grooves. This option has not been implemented at 
this time. 


7 



FUNCTIONING OF CODE 


This section describes in general terms how the code operates. Figure 3 shows a very simplified 
flow chart with only major subroutines. After the user calls up the code by typing LERCHP, subrou- 
tine INPUT is called. Data are entered from the keyboard as described in the following section. At 
the completion of data input, subroutine RUNCAS is called. This subroutine controls the flow of exe- 
cution. First it calls subroutine GUESS, in which an initial guess is made of the evaporator surface 
temperature at the beginning of the pipe T4E. The code uses all the heat input, environment tempera- 
ture, radiation parameter, and heat transfer coefficient information that has been provided in INPUT to 
make this guess. RUNCAS then calls subroutine START, which initializes other temperatures at the 
various radial locations at the beginning of the pipe. If the estimate of T4E is near the critical 
temperature of the working fluid at the outset, the user is returned to INPUT. 

Following this, RUNCAS begins an iteration to determine the operating condition of the heat 
pipe. Subroutine SUMEP is called. This computes the pressure, temperatures, and axial heat flux 
along the pipe by using a Runge-Kutta algorithm presented in reference 19. When the computation 
along the entire pipe length is completed, the axial heat flux at the downstream end of the pipe should 
be zero. If this is not the case, a correction to T4E is estimated and the iteration continues until the 
residual axial heat flux is less than 0.01% of the heat input. During each iteration the flow is 
monitored for the sonic condition. If this is encountered before reaching the downstream end of the 
last heat input or adiabatic section, T4E is raised to the minimum value needed to move this sonic 
point to the end of the section. The iteration is begun again at the beginning of the pipe. The flow is 
monitored for the viscous limit by comparing the vapor pressure computed at any point to the pressure 
at the evaporator end cap. If the ratio is less than 0.1, execution is terminated and a message is 
printed instructing the investigator to provide new input conditions. Alternatively, a viscous limit is 
deemed to be present if no temperature can be found at which the entire condenser length will reject 
heat. 


After convergence is reached, the data print subroutine is called. Data are printed out at the 
chosen axial positions along the pipe, and at the beginning and end of each pipe section. A new data 
input is then requested. 


Description of Data Input 

In the present PC version of the program, parameter input is by means of an interactive 
subroutine. The user is queried concerning pipe geometry, wick data, and choice of working fluid. 

This section demonstrates the operation of the code. Note that in this and following sections of this 
manual, statements written to the screen or to the printer from the code are written in upper case. 
Narrative comments about these statements are written in lower case. These comments only appear 
herein and are not printed by the code. The input data for an example case, without the descriptive 
material, appear in appendix C. 

The executable version of the code is file LERCHP.EXE, which stands for Lewis Research 
Center Heat Pipe. The suffix EXE means that the file is an executable file. The executable version of 
the code is started by typing LERCHP at the system operating level. The code will respond by 
typing the following statement on the screen and the printer: 
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LEWIS HEAT PIPE CODE, NASA LEWIS RESEARCH CENTER 

The code prints out input data to the printer periodically as the user goes through the input subroutine. 

IS THE HEAT PIPE OPERATING IN A “G” FIELD? TYPE “Y” 

The code asks the user if the heat pipe is operating in a G field. If the answer to this question is no, 
hit the Enter key or answer with N. The code will then assume that the heat pipe is operating in zero 
gravity. If the answer to this question is yes, the following message is printed: 

“G” FIELD STRENGTH, FRACTIONS OF “G” 0.00000000 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

The G fields for heat pipes operating on Earth, the Moon, and Mars are 1.0. 0.17, and 0.33, 
respectively. 

BEGIN INPUT OF HEAT PIPE GEOMETRY 

NUMBER OF RUNGE-KUTTA STEPS CHOSEN. 100 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

The number entered here will correspond to calculation points along the heat pipe axis. The default 
value is 100, but more Runge-Kutta steps may be introduced as needed to avoid large changes in 
conditions from one point to the next. 

LIMIT ON DATA POINTS PRINTED IS 100. 

SHOULD NOT EXCEED NUMBER OF RUNGE-KUTTA STEPS. 

NUMBER OF DATA POINTS PRINTED ALONG THE PIPE 50 
IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

The default value for print output of calculation results is 50 out of a possible 100 points. 

NUMBER OF PIPE SECTIONS TO BE CONSIDERED 2 
IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

This is the total number of heat pipe sections. Sections can be defined as regions of differing thermal 
boundary conditions or differing structure. The default value is 2. 

ARE PIPE DIMENSIONS AND WICKING SAME IN ALL SECTIONS? TYPE “Y” 

If yes is answered to this question, wick properties for all sections are assumed to be constant, and 
wick data need only be entered for the first section. The code will assign the same input data to all 
other sections. If the wick structure or containment envelop change from one axial location of the 
heat pipe to another, do not type Y but press the enter key. In this case the questions listed below will 
reappear for each section of the pipe. If input is repeated, Y must be typed in again if the wicking in 
the pipe sections is the same. 

If yes is answered to the preceding question, the following will appear on the screen: 
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WICKING AND DIAMETERS ALL SAME AS SECTION 1. 
OUTSIDE PIPE WALL DIAMETER, DO, m .00000000 
IF CHANGE, TYPE “Y”, THEN PRESS ENTER 


Input the outside diameter of the containment envelope in meters. 


INSIDE PIPE WALL DIAMETER, DI, m .00000000 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 


Input the inside diameter of the containment envelope in meters. 

WICK TYPE CODE 0 

1 = USER INPUT DATA 

2 = SCREEN WICK 

3 = SCREEN WITH ARTERIES 

4 = SINTERED METAL 

5 = SINTERED METAL WITH ARTERIES 

6 = GROOVES 

CHOICE OF WICK CODE 0 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 


There are six wick options in the code at this time. There is a section in the manual that describes the 
input data required for each case. If choice 6 is requested, the user is reminded that this choice is not 
implemented. In this case the user input option 1 was chosen and the data required are 

USER INPUT WICK PROPERTIES. 

WICK PERMEABILITY * AREA, AK, m**4 .00000000 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

This is the permeability of the wick to fluid flow in square meters multiplied by the flow area of the 
wick in square meters. 

VAPOR SPARE PERIMETER, PE, m .00000000 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

This is the perimeter at the liquid-vapor interface in meters. 


POROSITY OF WICK, PORE 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 




Input the volume of liquid in the wick divided by the total volume of the wick. 

EFFECTIVE WICK CONDUCTIVITY, KE, Wlm K .00000000 
IF CHANGE, TYPE “Y”, THEN PRESS ENTER 


This is the effective thermal conductivity of the wick filled with fluid, in the radial direction. 

EFFECTIVE VAPOR AREA, AV, m**2 .00000000 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 
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Input the effective cross-sectional area perpendicular to the axis. 


EFF. MIN CAPILLARY RADIUS, RCM, m 
IF CHANGE, TYPE “Y”, THEN PRESS ENTER 




This is the effective minimum radius of curvature of the meniscus, which is used in calculating the 
capillary limit. 


ENTRAINMENT DIMENSION, DCRITE, m 
IF CHANGE, TYPE “Y”, THEN PRESS ENTER 




Input the critical dimension used to calculate the entrainment limit. See reference 20. 

CRITICAL NUCLEATION RADIUS, RN, m .00000000 

DEFAULT VALUE = 1.0E-6 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

The minimum nucleation radius is a parameter used to calculate the boiling limit. This parameter is a 
function of the surface finish on the inside of the containment envelop. The type of wick structure 
also has some effect. The boiling limit equation and a description of this parameter can be found in 
reference 2 1 . 

TO REVIEW OR CORRECT THE INPUT, TYPE “Y” 

If a mistake is made when typing input into the code or if the user wants to review input for errors 
before a run, typing Y will allow the user to go back to the beginning of the input. The data values 
previously entered will be printed out on the screen next to the questions. 

Here a summary of the data input up to this point will be output by the printer. 

MATERIAL CHOICE CODE 

1 = STAINLESS STEEL 304 

2 = TITANIUM 

3 - IRON 

4 = NICKEL 

5 - ALUMINUM 

6 = COPPER 

7 = NIOBIUM-1 % ZIRCONIUM 

8 = TUNGSTEN 

9 = INCONEL 

10 = MOLYBDENUM 

11 = TANTALUM 

12 = BERYLLIUM 

TYPE OF PIPE MATERIAL 0 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

Enter the code from the preceding list that corresponds to the material of the heat pipe containment 
envelope. 
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TYPE OF WICK MATERIAL 0 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 


Enter the code from the preceding list that corresponds to the material of the heat pipe wick. 

FLUID TYPE CODE 

1 = LITHIUM 

2 = POTASSIUM 

3 = SODIUM 

4 = WATER 

5 = AMMONIA 

6 = METHANOL 

7 = MERCURY 

8 = DOWTHERM A 

9 = ACETONE 

10 = BENZENE 

11 = CESIUM 

12 = FREON-11 

13 = N-HEPTANE 

14 = SILVER 

15 = TOLUENE 

16 = M-XYLENE 

FLUID TYPE 0 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

Enter the code from this list that corresponds to the heat pipe working fluid. 


WETTING ANGLE, deg .00000000 
IF CHANGE, TYPE “Y”, THEN PRESS ENTER 


Enter the wetting angle, in degrees, of the fluid on the wick material. This value is typically assumed 
to be 0.0. In a heat pipe the wetting angle must be less than 90°. 

WEBER NUMBER, DEFAULT 2*PI, SEE INSTRUCTIONS 6.2832 
IF CHANGE, TYPE “Y”, THEN PRESS ENTER 


The Weber number is used to determine if the entrainment limit has been reached or exceeded. The 
Weber number at which entrainment occurs is not well understood in theory and may be different for 
different heat pipe wick structures. A value of 27T is set as the default. This was taken from reference 
20 . 

AVAILABLE CHOICES OF BOUNDARY CONDITIONS FOR SECTIONS 

1 HEAT INPUT OR OUTPUT SPECIFIED IN SECTION 

2 SURFACE TEMPERATURE SPECIFIED IN SECTION 

3 ENVIRONMENT TEMPERATURE GIVEN IN SECTION 

CHOICE FOR SECTION # 1 >0 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 


12 


Choose the boundary condition code that corresponds to the case of interest. This message will be 
printed out once for each section of the pipe. Option 1 is used when specified heat inputs and outputs 
are known. Option 2 can be used when surface temperatures and heat transfer data to or from the 
environment are known, but option 2 currently is not operational. Option 3 is used when the 
environment temperature and the heat transfer data between the environment and the pipe are 
specified. 

The following input routine is for option 1. 

HEAT IN OR OUT ALONG SECTION IS CONSTANT 
IF CHANGE, TYPE “Y” 

If the heat input to a section, section 1 in this case, is constant along that section, do not type Y, just 
hit the enter key. Doing this causes the code to assume a constant heat input along the section. The 
code will then prompt the user for the following: 

AXIAL LOCATION OF SECTION BEGINNING, m .00000000 
IF CHANGE, TYPE “Y” 

Zero is normally taken as the upstream end of the evaporator, but it is possible to designate a nonzero 
value of the axial coordinate as the beginning of the pipe. 

AXIAL LOCATION OF SECTION END, m .00000000 
IF CHANGE, TYPE “Y” 

This position is measured with respect to the chosen zero of the coordinate system. 

HEAT INTO OR OUT OF SECTION, W .00000000 

IF CHANGE, TYPE “Y M 

For sections with constant heat input or output, heat is entered in watts. If heat input is not constant 
along the section of interest, the user should type Y to the query IF CHANGE, TYPE “Y” following 
HEAT IN OR OUT ALONG SECTION IS CONSTANT, thusly: 

HEAT IN OR OUT ALONG SECTION IS CONSTANT 
IF CHANGE, TYPE “Y” 

Type Y here and the code will respond with 

YOU HAVE CHANGED THE CHOICE TO VARIABLE 

NUMBER OF SECTION POINTS CONSIDERED > 0 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

The user must enter the number of points along the section where the heat input or output rate is 
known. After this is done, the code will then ask the user the location and heat transfer rate at each 
point. 

POINT# 1 
Z( 1) .00000000 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 
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Input the location where the heat input rate is known. This position is measured with respect to the 
upstream end of the evaporator. 

DQI( 1) .00000000 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

Input the heat load in watts per meter of pipe length. It is important to note that an axially varying 
heat load is entered in watts per meter. However, for sections with constant heat input or output, heat 
is entered in watts. 

If option 3 is chosen, the code will respond 

ENVIRONMENT TEMPERATURE, K .0000000 
IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

The environment temperature should be entered here. For a radiation heater this would be the heater 
surface temperature. For a heat pipe that is cooled by radiation heat transfer to space, this temperature 
would be the effective environment temperature of deep space. 

AXIAL LOCATION OF SECTION BEGINNING, m .00000000 
IF CHANGE, TYPE “Y” 

Zero is normally taken as the upstream end of the evaporator, but it is possible to designate a nonzero 
value of the axial coordinate as the beginning of the pipe. 

AXIAL LOCATION OF SECTION END, m .00000000 
IF CHANGE, TYPE “Y” 

This position is measured with respect to the chosen zero of the coordinate system. 

SPECIFIC AREA PER UNIT LENGTH, SA .00000000 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

This is the area in square meters per meter of axial distance from which the heat pipe either absorbs or 
radiates heat. If a heat pipe has fins attached to the condenser section, the specific radiating area 
would be the area of the fins plus the area of the pipe for 1 m of length. 

RADIATION FACTOR, RADK .00000000 

IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

The radiation factor RADK is used to account for the effective emissivity of the pipe and fin surfaces. 
This parameter is also used to account for overall view factors and fin efficiency between the heat pipe 
and fin surfaces. 

CONVECTION COEFFICIENT, HC, Wfm**2*K .00000000 
IF CHANGE, TYPE “Y”, THEN PRESS ENTER 

This is the value of the convection coefficient between the heat pipe wall and the environment 
temperature input. 
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After boundary conditions have been entered for all the sections of the pipe, the following message 
will appear: 

TO REVIEW OR CORRECT MATERIALS, FLUID, BOUNDARY CONDITIONS, TYPE “Y” 

If a mistake is made when typing input in the code or if the user wants to review input for errors 
before a run, typing “Y” here will allow the user to go back to the point at which material choices are 
inserted. The data values previously entered will be printed out on the screen next to the questions. 

After this point a summary of the data input up to this point will be output by the printer. After the 
input has been printed, the code will respond with 


COMPUTING > 


Description of Data Output 

After a solution has been obtained, the operator is queried: 

DO YOU WISH TO OMIT PRINTOUT OF DATA ALONG THE PIPE? TYPE Y 

Data will not be printed out along the pipe. Rather, a summary of the upstream evaporator and 
downstream condenser temperatures will be printed. Otherwise the printer will send the message 

BEGIN PRINTOUT OF RESULTS 

After this message, data will be printed at the user-specified points. Variable names and the 
corresponding numeric values are printed. Dimensions are omitted. A table describing the variable 
names and dimensions is given in appendix A. Also see appendixes D and E for example output 
cases. 

Output will be generated at the printer for each run. It includes the following: 

(1) All user inputs 

(2) Calculated heat transport 

(3) Liquid pressures 

(4) Vapor pressures 

(5) Interpretive messages 

Values for the following parameters are given for each printout point in a section: 


(1) Axial location 

(2) Liquid and vapor pressures 

(3) Liquid-vapor interface, outer surface, and mean vapor temperatures 

(4) Mean vapor velocity 

(5) Mean Mach number 

(6) Velocity profile factor 

(7) Radial heat input or output per unit length 

(8) Axial vapor heat flux 
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APPENDIX A 


A 

AK 

AV 

DELTCR 

DCRITE 

DI 

DO 

DQDZ 

HC 

KE 

MACHM 

NA 

P 

PE 

PL 

PORE 

Q 

RADK 

RCM 

RN 

RSS 


SYMBOLS USED IN INPUT AND OUTPUT 
Busse velocity profile parameter 

product of wick cross-sectional area and permeability, m 4 
cross-sectional area of vapor passage, m 2 

temperature drop through evaporator wick for initiation of boiling, K 
entrainment critical dimension, m 
inner diameter of pipe wall, m 
outer diameter of pipe wall, m 

radial heat input into or out of pipe per unit length, W/m 
heat transfer coefficient, W/m 2 K 

effective thermal conductivity of liquid-saturated wick, W/m K 

mean Mach number at axial position Z computed by using the mean velocity VM, with 
sound speed computed from mean temperature TM 

number of arteries 

vapor pressure at axial position Z, N/m 2 

perimeter of vapor passage cross section, m 

liquid pressure in wick at position Z, N/m 2 

porosity of wick, volume of void in wick per wick volume 

axial vapor heat transport at position Z, W 

nondimensional factor including emissivity, view factor, and fin efficiency 
effective minimum capillary radius in wick, m 
minimum nucleation radius, m 
radius of particles in sintered wick, m 
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SA 

TBOIL 


TM 

TWI 

T1Z 

T2Z 

T3Z 

T4Z 

VM 

Z 


specific area, m 2 /m 

temperature drop through evaporator wick, which will result in boiling if DELTCR is 
exceeded, K 

mean temperature at axial position Z, K 
wick thickness, m 

temperature of vapor at vapor-liquid interface at axial position Z, K 
temperature of liquid at vapor-liquid interface at axial position Z, K 
temperature of inner surface of pipe at axial position Z, K 
temperature of the outer surface of pipe at axial position Z, K 
mean velocity at axial position Z, m/s 
axial distance from beginning of pipe, m 


The fluid codes are as follows: 


1 

lithium 

9 

acetone 

2 

potassium 

10 

benzene 

3 

sodium 

11 

cesium 

4 

water 

12 

Freon- 1 1 

5 

ammonia 

13 

n-heptane 

6 

methanol 

14 

silver 

7 

mercury 

15 

toluene 

8 

Dowtherm A 

16 

m -xylene 


The material codes are as follows: 

1 stainless steel 304 

2 titanium 

3 iron 

4 nickel 

5 aluminum 

6 copper 

7 niobium-1 % zirconium 

8 tungsten 

9 Inconel 

10 molybdenum 

1 1 tantalum 

12 beryllium 
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The wick codes are as follows: 

1 user input data 

2 screen wick 

3 screen wick with arteries 

4 sintered metal 

5 sintered metal with arteries 

6 grooves 



APPENDIX B 


VAPOR FLOW ALGORITHM* 

Analysis 

In the analysis to follow, the appropriate features of Busse’s (ref, 16) method will be 
introduced, and the development will proceed from there. In addition, some major assumptions for the 
heat pipe sections being considered should be indicated at the outset: Flow is laminar. The ideal-gas 
law applies to components of the vapor. The equations allow for local chemical equilibrium among 
species in the vapor. Phase equilibrium occurs only at the liquid-vapor interface, so that liquid 
droplets are absent in the vapor. Following Busse (ref. 16), no tangential velocity component exists at 
the wall in mass addition or removal, the radial pressure gradient may be neglected, and the velocity 
components are zero at the closed ends of the heat pipe. Also, his method of imposing a negligible 
radial pressure gradient by equating the mean axial pressure gradient over the cross section to the axial 
pressure gradient along the centerline was used. 

Busse (ref. 16) began with the axial axisymmetric Navier-Stokes equation 



where p is pressure, r is radial distance, u is radial velocity, v is axial velocity, z is axial distance, r| 
is dynamic viscosity, and p is density. Busse integrated the mass conservation equation with respect 
to «, so that this term can be eliminated. This is also done in the present analysis except that density 
is variable in z and n 


u = 



( 2 ) 


When equation (2) is combined with equation (1), the following results: 


dp _ 1 
dz r 




„„ dv n d dv 

- pv — + _ — r 

dz r dr dr 


( 3 ) 


Here the axial variation of u has been neglected in the viscous term. The assumption is made that the 
radial variation of density is much less than the radial variation of axial velocity because the latter 
vanishes at the liquid-vapor interface. On the basis of this assumption, a presently undefined mean 
density p m and its axial derivative are extracted from the integrals, yielding 


dp_ 

dz 


dv 

Jri K 


r 


dv 

dz 


r' dr ' + 


dz 


r 


vr' dr' 





q d r dv 
r dr dr 


( 4 ) 


‘This appendix is taken from reference 17, except for the last two paragraphs, which are new. 
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The crux of Busse’s attack on the incompressible case is the assumption of an approximate velocity 
profile that satisfies all the boundary conditions. It incorporates a parameter A(z) that may be constant 
or may vary with position, depending on mass influx boundary conditions. The same profile is 
employed in this analysis: 


f 

v — 2v 1 

m 

\ 



I 


1 + A(z) . 

\ 




(5) 


Here v m is the mean axial velocity and R is the vapor space radius. When A(z) — 0, equation (5) 
reduces to the Poiseuille profile. Integrating equation (4) over the cross section of the pipe and 
manipulating give 



( 6 ) 


The mean density is assumed to be related to the pressure and a mean temperature T m by the ideal- 
gas law, which yields upon differentiation 


d ln Pm 
dz 


= (GRP) fi lg 

dz 


- (GRT) 


d In T 

m 

dz 


(7) 


where 


5 In R 

GRP - 1 - 1 and 

0 In p 


GRT = 1 


0 In R 

+ £ 

0 In T 

m 


and where R g , the gas constant, may depend on pressure and temperature if chemical reactions are 
considered. 


Continuity requires that the vapor mass flux m(z) at any cross section equal the product of mean 
density p m , the vapor cross-sectional area A v , and the mean axial velocity v m (z). From this relationship 
and equation (7) the axial derivative of v m is found to be 

d ln v A In m d In A A In n d \n T 

» = a in m - l - (GRP) a m P + (GRT) * ( 8 ) 

dz dz dz dz dz 


Equations (6) to (8) are combined and the indicated integrations are performed. Dropping the 
distinction between the mean cross-sectional pressure and the local pressure by Busse’s approximation 
gives 


1 


4 Pm V « 

3p 


(GB) (GRP) 


d ln p 
dz 


4p v 

~ Ml W) 


3p 


(GB)(GRT) 


d In T 

m 

dz 


4p 

m m 

~^P~ 


(HB) 


dA 

dz 


00 

fl 

8p vl 

- ” m (OR) 

( 

d In m 

\ 

d ln A 

__ V 

R 2 p 

l 3 J 

3/> 

dz 

dz J 


( 9 ) 
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where 


GB = 

fi-4 

2A 2 ^ 

+ 

and 

HB = - 

( 1 4 AX 


l 6 

45 

J 



[6 " 45j 


The variation of A v with z is understood to be small. 


The mean temperature T m is evaluated by means of an energy integral. Figure 4 shows the control 
volume through which energy is transported. The energy flows into the boundary are equated with the 
energy flows out, internal energy and pressure are eliminated by using enthalpy, and the limit Az — > 0 
is taken, to give 


f 


dz 


pv 


\ + — 




r dr = <Fr dr + qjt - pRu R 


( 10 ) 


where h is the enthalpy of the mixture per unit mass, q R is the energy flux entering the vapor space 
radially, and <E> is the viscous dissipation per unit volume. The relation for u„, the radial velocity at 
the vapor-liquid interface, is given by the mass conservation equation in the form 


«* - 



( 11 ) 


The energy influx q R is specified as -p R u R (e R + ify 2), where e R is internal energy. The result of 
combining equations (10) and (11) is 



v 

T 


H. 


-jL (pv) r dr + 


dh . 
pv — r dr + 
dz 


f 


pv* 


dv 

dz 


r dr 



< 12 ) 


The mean value of both the enthalpy and the density can be separately extracted from the integrals 
if the assumption is again made that the variation of both variables with radius is much less than that 
of axial velocity. When these mean values have been extracted and terms common to all expressions 
cancelled, the energy equation becomes 


K 




X 


*dv 3 

dz 


r dr - 


2 n 
~R* Jo 


<Fr dr 


where 


(13) 


If chemical equilibrium in the vapor must be considered, enthalpy h m becomes a function of both 
T m and p. When its derivative is inserted into equation (13) with equations (7) and (9) and the 
integrations are performed, the result is 
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( \ 
dh i 


d In p 


- FV (GRP) 
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d In p + 

dh m 

dz 

d In T 

m Jp 


+ F 3 v m (GRT) 


\d In T dA 

2 + (HAF) _ 
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P 
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dz 


dz 
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where 


F 3 - 8 


HAF = - 


< 1 

V' 
2 ( 


A + A 2 2A 3 
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An additional equation required to find the three dependent variables p. A, and T m comes from 
Busse’s (ref. 16) approximation that the axial gradient of the mean pressure equals the centerline 
pressure gradient. The latter, from equation (1) for r = 0 using equations (4) and (5), is 


A 2 

1 - PmVm (1 - A/3) 2 (GRP) 

d In p + 

A 2 

PmVm (1 - A/3) 2 (GRT) 
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dz 
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m 

dz 
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dz 
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(15) 


The mass flux is related to the heat flux into the pipe by 


1 dm 


d( lR 


dz 2nR dz 


( 


K, + 




1 dm 


27T/?p R dz 


J J 


(16) 


where the kinetic energy of the injected vapor has been incorporated, and h vl is the heat of evaporation. 
Equations (9), (14), and (15), with equation (16) and appropriate boundary conditions, enable a solu- 
tion for p. A, and T m . 

At extremely low Mach numbers the flow approaches the incompressible condition. The solution 
algorithm defined by the foregoing equations breaks down because of roundoff error in the computer. 
A simple fix can be devised for the very low peak Mach number regime by returning to the original 
equations of Busse (ref. 16) derived for incompressible conditions. 
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Busse’s equations (18) and (20) can be combined to obtain an expression for the change of mean 
pressure with distance in the following form: 


dp 

dz 


pv„(180 - 48A + 4A 2 ) 


d In v 8nv 

m + * m 

dz R 2 

75 - 22A 


75 + 58A - 





1A 


(17) 


In the equation the mean velocity derivative is eliminated by the substitution 

d In m u d ln v m + d In A y 
dz dz dz 


The change of velocity profile parameter can be expressed as 


dA 

dz 


(90 - 150A + 22A 2 ) _ - 1 040,1 A 

dz R 2 pv 

r m 

' 75 - 22A 


(18) 


In the code equations (17) and (18) are solved for each step as though the flow were incompressible. 
However, the density is updated by means of the gas law. The saturated vapor interface temperature is 
employed in the gas law. This procedure was used when the highest Mach number at any point in the 
pipe fell below 0.01. 
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APPENDIX C 


USER-INTERACTIVE DATA ENTRY FOR FIRST EXAMPLE CASE 

The example case given here demonstrates the flexibility of this code. The heat pipe is made of 
stainless steel and has sodium as a working fluid. The heat pipe has four different sections: two 
evaporators, one adiabatic section, and one condenser section. The heat input profiles are known for 
both evaporators and are shown graphically in figure 5. The heat input is constant for the upstream 
evaporator but variable for the other evaporator. The heat input is known at three points along the 
variable evaporator section. The code will solve for the temperature distribution along the heat pipe 
and the heat load distribution in the condenser section. 

The input procedure for the example given here is exactly as it appears on the computer screen. 

LEWIS HEAT PIPE CODE, NASA LEWIS RESEARCH CENTER 

BEGIN INPUT OF HEAT PIPE GEOMETRY 

IS THE PIPE IN A "G" FIELD? TYPE "Y" 


INITIAL NUMBER OF RUNGE-KUTTA STEPS CHOSEN 100 
IF CHANGE, TYPE "Y", THEN PRESS ENTER 


LIMIT ON DATA POINTS PRINTED ALONG PIPE IS 100. 

SHOULD NOT EXCEED NUMBER OF RUNGE-KUTTA STEPS. 

NUMBER OF DATA POINTS PRINTED ALONG THE PIPE 50 
IF CHANGE, TYPE ,, Y H , THEN PRESS ENTER 

Y 

ENTER NEW INTEGER VALUE 
40 

NUMBER OF PIPE SECTIONS TO BE CONSIDERED 2 

IF CHANGE, TYPE "Y" , THEN PRESS ENTER 

Y 

ENTER NEW INTEGER VALUE 
4 

ARE PIPE DIMENSIONS AND WICKING SAME IN ALL SECTIONS? TYPE "Y” 

Y 

OUTSIDE PIPE WALL DIAMETER, DO, m .00000000 
IF CHANGE, TYPE "Y M , THEN PRESS ENTER 

Y 

ENTER NEW REAL VALUE 
.0191 

INSIDE PIPE WALL DIAMETER, DI, m .00000000 
IF CHANGE, TYPE M Y” , THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
.0175 

WICK TYPE CODE 

1 = USER INPUT DATA 

2 = SCREEN WICK 

3 = SCREEN WITH ARTERIES 

4 * SINTERED METAL 

5 = SINTERED METAL WITH ARTERIES 

6 = GROOVE 

CHOICE OF WICK CODE 0 

IF CHANGE, TYPE M Y", THEN PRESS ENTER 

Y 

ENTER NEW INTEGER VALUE 

1 

WICK PERMEABILITY* AREA, AK, m**4 .00000000 

IF CHANGE, TYPE "Y'\ THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
4 . 87e-14 
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VAPOR SPACE PERIMETER, PE, m .00000000 

IF CHANGE, TYPE "Y", THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
.0449 

POROSITY OF WICK, PORE .00000000 

IF CHANGE, TYPE "Y", THEN PRESS ENTER 
Y 

ENTER NEW REAL VALUE 
.65 

EFF. WICK CONDUCTIVITY, KE, W/m K .00000000 
IF CHANGE, TYPE "Y" , THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
46 

EFFECTIVE VAPOR AREA, AV, m**2 .00000000 

IF CHANGE, TYPE "Y M , THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
1. 606e-4 

EFF. MIN CAPILLARY RADIUS, RCM, m .00000000 
IF CHANGE, TYPE "Y M , THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
3 . 81e-5 

ENTRAINMENT DIMENSION, DCRITE, M .00000000 

IF CHANGE, TYPE H Y", THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
3 . 8 le-5 

UNLESS WICK IS UNIFORM, 

PROVIDE WICK THICKNESS FOR BOILING 
LIMIT CALCULATION CALC, TWI , m .00000000 

IF CHANGE, TYPE M Y M , THEN PRESS ENTER 

NUCLEATION RADIUS, RN, TO .00000000 

DEFAULT VALUE - 1.0E-6 

IF CHANGE, TYPE "Y", THEN PRESS ENTER 

TO REVIEW OR CORRECT THE INPUT, TYPE M Y M 


WICKING , X-SECTION DIMENSIONS SAME AS SECTION 1 

MATERIAL CHOICE CODE 

1 = STAINLESS STEEL 304 

2 = TITANIUM 

3 = IRON 

4 - NICKEL 

5 = ALUMINUM 

6 = COPPER 

7 = NIOBIUM- 1 % ZIRCONIUM 

8 = TUNGSTEN 

9 = INCONEL 

10 = MOLYBDENUM 

11 = TANTALUM 

12 = BERYLLIUM 

TYPE OF PIPE MATERIAL 0 

IF CHANGE, TYPE "Y H , THEN PRESS ENTER 

y 

ENTER NEW INTEGER VALUE 

1 

TYPE OF WICK MATERIAL 0 

IF CHANGE, TYPE "Y", THEN PRESS ENTER 

y 

ENTER NEW INTEGER VALUE 

1 
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FLUID TYPE CODE 0 

1 = LITHIUM 

2 = POTASSIUM 

3 - SODIUM 

4 = WATER 

5 = AMMONIA 

6 = METHANOL 

7 * MERCURY 

8 = DOWTHERM A 

9 = ACETONE 

10 = BENZENE 

11 = CESIUM 

12 - FREON-11 

13 = N-HEPTANE 

14 = SILVER 

15 = TOLUENE 

16 = M-XYLENE 

FLUID TYPE 0 

IF CHANGE , TYPE "Y" , THEN PRESS ENTER 

Y 

ENTER NEW INTEGER VALUE 
3 

WETTING ANGLE, DEG .00000000 

IF CHANGE, TYPE "Y", THEN PRESS ENTER 

WEBER CONSTANT, DEFAULT 2*PI,SEE INSTRUCTIONS 6.2831800 
IF CHANGE, TYPE M Y M , THEN PRESS ENTER 

AVAILABLE CHOICES OF BOUNDARY CONDITIONS FOR SECTIONS 

1 HEAT INPUT OR OUTPUT SPECIFIED IN SECTION 

2 SURFACE TEMPERATURE SPECIFIED IN SECTION 

3 ENVIRONMENT TEMPERATURE GIVEN IN SECTION 

CHOICE FOR SECTION # 1 >0 

IF CHANGE, TYPE "Y" , THEN PRESS ENTER 

y 

ENTER NEW INTEGER VALUE 

1 

HEAT IN OR OUT ALONG SECTION IS CONSTANT 
IF CHANGE, TYPE M Y M 

AXIAL LOCATION OF PIPE BEGINNING, M .00000000 

IF CHANGE, TYPE M Y H , THEN PRESS ENTER 

AXIAL LOCATION OF SECTION END, M .00000000 

IF CHANGE, TYPE "Y", THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 

.2 

HEAT INTO OR OUT OF SECTION , WATTS .00000000 

IF CHANGE, TYPE M Y" , THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
2000 

AVAILABLE CHOICES OF BOUNDARY CONDITIONS FOR SECTIONS 

1 HEAT INPUT OR OUTPUT SPECIFIED IN SECTION 

2 SURFACE TEMPERATURE SPECIFIED IN SECTION 

3 ENVIRONMENT TEMPERATURE GIVEN IN SECTION 

CHOICE FOR SECTION # 2 >0 

IF CHANGE, TYPE H Y" , THEN PRESS ENTER 

y 

ENTER NEW INTEGER VALUE 

1 

HEAT IN OR OUT ALONG SECTION IS CONSTANT 
IF CHANGE, TYPE "Y" 

AXIAL LOCATION OF SECTION END, M .00000000 

IF CHANGE, TYPE "Y M , THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
.3 

HEAT INTO OR OUT OF SECTION , WATTS .00000000 

IF CHANGE, TYPE "Y" , THEN PRESS ENTER 
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AVAILABLE CHOICES OF BOUNDARY CONDITIONS FOR SECTIONS 

1 HEAT INPUT OR OUTPUT SPECIFIED IN SECTION 

2 SURFACE TEMPERATURE SPECIFIED IN SECTION 

3 ENVIRONMENT TEMPERATURE GIVEN IN SECTION 

CHOICE FOR SECTION # 3 >0 

IF CHANGE, TYPE "Y" , THEN PRESS ENTER 

y 

ENTER NEW INTEGER VALUE 

1 

HEAT IN OR OUT ALONG SECTION IS CONSTANT 
IF CHANGE, TYPE "Y" 

Y 

YOU HAVE CHANGED THE CHOICE TO VARIABLE. 

NUMBER OF SECTION POINTS CONSIDERED > 0 

IF CHANGE, TYPE M Y M , THEN PRESS ENTER 

y 

ENTER NEW INTEGER VALUE 
3 

POINT # 1 

Z( 1) .30000000 

IF CHANGE, TYPE "Y", THEN PRESS ENTER 

DQI ( 1) .00000000 

IF CHANGE, TYPE "Y", THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 

10000 
POINT # 2 

Z( 2) .00000000 

IF CHANGE, TYPE M Y", THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 

.4 

DQI ( 2) .00000000 

IF CHANGE, TYPE "Y M , THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
11000 
POINT # 3 

Z( 3) .00000000 

IF CHANGE, TYPE ”Y", THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
.5 

DQI ( 3) .00000000 

IF CHANGE, TYPE "Y M , THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
12500 

AVAILABLE CHOICES OF BOUNDARY CONDITIONS FOR SECTIONS 

1 HEAT INPUT OR OUTPUT SPECIFIED IN SECTION 

2 SURFACE TEMPERATURE SPECIFIED IN SECTION 

3 ENVIRONMENT TEMPERATURE GIVEN IN SECTION 

CHOICE FOR SECTION # 4 >0 

IF CHANGE, TYPE "Y", THEN PRESS ENTER 

y 

ENTER NEW INTEGER VALUE 
3 

ENVIRONMENT TEMPERATURE, K .00000000 

IF CHANGE, TYPE "Y", THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
273 

DIMENSIONS, PROPERTIES IN SECTION ARE CONSTANT 
IF CHANGE, TYPE "Y" 

AXIAL LOCATION OF SECTION END, M .00000000 

IF CHANGE, TYPE M Y" , THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 

1.2 

SPECIFIC AREA PER UNIT LENGTH, SA .00000000 
IF CHANGE, TYPE "Y", THEN PRESS ENTER 

y 

ENTER NEW REAL VALUE 
.15 
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RADIATION FACTOR RADK .00000000 

IF CHANGE, TYPE "Y" , THEN PRESS ENTER 
Y 

ENTER NEW REAL VALUE 

.8 

CONVECTION COEFFICIENT HC, W/m**2 K .00000000 
IF CHANGE, TYPE "Y M , THEN PRESS ENTER 

TO REVIEW OR CORRECT MATERIALS, FLUIDS, BOUNDARY CONDITIONS, TYPE "Y n 
COMPUTING > 

UPSTREAM SURFACE T4E= 986.20090 
COMPLETED PASS 1 
COMPLETED PASS 1 

UPSTREAM SURFACE T4E= 993.54890 
COMPLETED PASS 2 
COMPLETED PASS 2 

UPSTREAM SURFACE T4E= 993.40250 

DO YOU WISH TO OMIT PRINTOUT OF DATA ALONG THE PIPE? TYPE Y 
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APPENDIX D 


OUTPUT DATA FOR EXAMPLE CASE 

The output data for the example case described in appendix C are given here. They are presented 
exactly as they appear by running the code. 


LEWIS HEAT PIPE CODE, NASA LEWIS RESEARCH CENTER 
BEGIN LISTING OF INPUT DATA 

"G" FIELD STRENGTH, FRACTIONS OF "G” .00000000 

INITIAL # OF RUNGE-KUTTA STEPS CHOSEN 100 
NUMBER OF DATA POINTS PRINTED ALONG PIPE 40 
NUMBER OF SECTIONS IN THE PIPE 4 

USER INPUT WICK PROPERTIES 

WICKING, X-SECTION DIMENSIONS SAME AS SECTION 1 

SECTION NUMBER 1 
USER INPUT DATA 

AK= . 48700000E-13 PE= . 44900000E-01 PORE= .65000000 

KE= 46.000000 AV= . 16060000E-03 DCRITE= . 38100000E-04 

RCM= . 38100000E-04 DO= . 19 100000E-01 DI= . 17500000E-01 

RSS= .00000000 RN= . 10000000E-05 TWI= .00000000 NA= 0 

TYPE OF PIPE MATERIAL CHOSEN 
1 = STAINLESS STEEL 304 
TYPE OF WICK MATERIAL CHOSEN 
1 = STAINLESS STEEL 304 

WEBER NUMBER 6.2831800 

WORKING FLUID 
SODIUM 

SECTION NUMBER 1 HEAT INPUT SPECIFIED IN THIS SECTION 

1=1 ZI= .00000000 1=2 Z 1= .20000000 QSEC= 2000.0000 

SECTION NUMBER 2 HEAT INPUT SPECIFIED IN THIS SECTION 

1=1 ZI= .20000000 1=2 ZI= .30000000 QSEC= .00000000 

SECTION NUMBER 3 HEAT INPUT SPECIFIED IN THIS SECTION 


# OF POINTS IN SECTION ( 3)= 3 


1= 1 

Z = 

.30000000 

DQI= 

10000.000 

1= 2 

Z = 

.40000000 

DQI= 

11000.000 

1= 3 

z= 

. 50000000 

DQI= 

12500.000 


SECTION NUMBER 4 ENVIRONMENT TEMPERATURE SPECIFIED IN THIS SECTION 
ENVIRONMENT TEMPERATURE , K 273.00000 
ZI ( 1) .50000000 ZI ( 2 ) 1.2000000 

SA= .15000000 RADK= .80000000 HC= .00000000 

UPSTREAM SURFACE T4E= 986.20090 
CALL SUMEP 
COMPLETED PASS 1 

UPSTREAM SURFACE T4E= 993.54890 
CALL SUMEP 
COMPLETED PASS 2 

UPSTREAM SURFACE T4E= 993.40250 
CALL SUMEP 

BEGIN PRINTOUT OF RESULTS 
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BEGIN SECTION # 1 
MENISCUS PRESS. DIFFERENCE 

VAPOR - LIQUID, N/SQ M * 3403.4530 

CAPILLARY LIMIT, N/SQ M * 6864.8180 


z* 

.00000000 


p* 

15478.560 

PL* 

12075.110 

T1Z* 

980.49690 


T4Z* 

993.40250 

TM= 

980.49690 

VM= 

.00000000 


MACHM= 

.00000000 

A* 

.54389510 

DQDZ= 

10000.000 


Q= 

.00000000 



Z= 

. 30000000E- 

01 

P= 

15473.500 

PL= 

12079.590 

T1Z* 

980.47120 


T4Z* 

993.37680 

TM= 

980.45940 

VM* 

9.7012210 


MACHM* 

. 15445590E-01 

A= 

.54396840 

DQDZ= 

10000.000 


Q= 

300.00000 



Z* 

. 60000000E- 

■01 

P* 

15457.290 

PL* 

12093.930 

T1Z— 

980.38870 


T4Z= 

993.29430 

TM* 

980.32920 

VM= 

19.419230 


MACHM* 

. 30920680E-01 

A* 

.54494620 

DQDZ= 

10000.000 


Q= 

600.00000 



Z- 

. 90000000E- 

■01 

P* 

15427.840 

PL* 

12119.930 

T1Z= 

980.23850 


T4Z= 

993.14410 

TM= 

980.09250 

VM= 

29.176580 


MACHM* 

. 46464980E-01 

A= 

.54616770 

DQDZ- 

10000.000 


Q= 

900.00010 



Z = 

.12000000 


P* 

15386.060 

PL* 

12156.700 

T1Z= 

980.02490 


T4 Z~ 

992.93050 

TM= 

979.75580 

VM* 

38.990310 


MACHM* 

. 62108500E-01 

A* 

.54788410 

DQDZ- 

10000.000 


Q* 

1200.0000 



Z= 

. 15000000 


P= 

15330.770 

PL* 

12205.140 

T1Z— 

979.74160 


T4Z= 

992.64730 

TM= 

979.30830 

VM* 

48.887650 


MACHM= 

. 778 9 923 0E- 01 

A* 

.55014200 

DQDZ = 

10000.000 


Q- 

■1500.0000 



Z= 

. 18000000 


P= 

15262.800 

PL* 

12264.360 

T1Z = 

979.39230 


T4Z* 

992.29800 

TM= 

978.75700 

VM= 

58.884090 


MACHM= 

. 93864520E-01 

A* 

.55291210 

DQDZ= 

10000.000 


Q“ 

1800.0000 



Z= 

.20000000 


P= 

15209.590 

PL* 

12310.440 

T1Z = 

979.11790 


T4Z* 

992.02350 

TM* 

978.32260 

VM= 

65.619220 


MACHM* 

.10463320 

A* 

. 55507700 

DQDZ* 

10000.000 


Q- 

1999.9810 



BEGIN 

SECTION # 2 






MENISCUS PRESS. DIFFERENCE 




VAPOR - LI QUID, N/SQ M = 

2899.1590 



CAPILLARY LIMIT, N/SQ 

M = 

6871.4060 



Z= 

.20000000 


P= 

15209.590 

PL= 

12310.440 

T1Z = 

979.11790 


T4 Z= 

979.11790 

TM* 

978 . 32260 

VM= 

65.619220 


MACHM* 

.10463320 

A* 

.55507700 

DQDZ= 

.00000000 


Q= 

1999.9810 



Z= 

.21000000 


P= 

15206.260 

PL* 

12335.370 

T1 Z= 

979.10060 


T4Z* 

979.10060 

TM* 

978.31480 

VM* 

65.633900 


MACHM* 

.10465630 

A= 

.54332270 

DQDZ= 

.00000000 


Q~ 

1999.9810 



Z= 

.24000000 


P= 

15196.370 

PL* 

12410.190 

T1Z= 

979.04960 


T4Z* 

979.04960 

TM* 

978.29270 

VM** 

65.677570 


MACHM= 

. 10472470 

A* 

.50979230 

DQDZ= 

.00000000 


Q“ 

1999.9810 



Z= 

.27000000 


P- 

15186.670 

PL* 

12485.010 

T1Z** 

978.99930 


T4Z- 

978.99930 

TM* 

978.27170 

VM* 

65.720570 


MACHM* 

.10479190 

A* 

.47874090 

DQDZ- 

.00000000 



1999.9810 



Z= 

.30000000 


P= 

15177.170 

PL* 

12559.830 

T1Z= 

978.95030 


T4Z* 

978.95030 

TM= 

978.25260 

VM* 

65.762900 


MACHM* 

.10485800 

A= 

.44999400 

DQDZ* 

.00000000 


Q= 

1999.9810 
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BEGIN SECTION # 3 
MENISCUS PRESS. DIFFERENCE 
VAPOR - LIQUID, N/SQ M = 
CAPILLARY LIMIT, N/SQ M = 


2617.3410 

6871.4060 


z= 

.30000000 

P= 

15177.170 

PL= 

12559.830 

T1Z= 

978.95030 

T4 Z — 

991.86660 

TM= 

978.25260 

VM= 

65.762900 

MACHM= 

.10485800 

A= 

.44999400 

DQDZ = 

10000.000 

Q= 

1999.9810 



Z— 

. 33000000 

P= 

15086.210 

PL= 

12638.270 

T1Z- 

978.47800 

T4Z= 

991.52550 

TM= 

977.48750 

VM= 

76.042410 

MACHM= 

.12131640 

A= 

. 48685010 

DQDZ= 

10101.520 

Q= 

2300.9140 



Z= 

.36000000 

P= 

14979.850 

PL= 

12727.610 

T1Z— 

977.92240 

T4Z= 

991.32700 

TM= 

976.59390 

VM= 

86.671330 

MACHM= 

.13836290 

A= 

. 51268230 

DQDZ= 

10378.000 

Q= 

2607.0600 



Z- 

. 39000000 

P= 

14853.530 

PL= 

12829 . 080 

T1Z= 

977.25820 

T4Z— 

991.23760 

TM= 

975.52530 

VM= 

97.914080 

MACHM= 

. 15643310 

A= 

.53274700 

DQDZ= 

10822.960 

Q= 

2924.2500 



Z— 

.42000000 

P= 

14704.910 

PL= 

12942 . 190 

T1Z= 

976.47110 

T4Z= 

991.19060 

TM= 

974.25710 

VM= 

109.94350 

MACHM= 

.17581380 

A= 

.55040510 

DQDZ= 

11396.000 

Q= 

3256.0190 



Z = 

.45000000 

P= 

14527.450 

PL= 

13068.600 

T1Z= 

975.52310 

T4 Z— 

990.99820 

TM= 

972.72350 

VM= 

123.02390 

MACHM- 

.19695300 

A= 

.56707200 

DQDZ = 

11981.000 

Q= 

3606.3140 



Z= 

.48000000 

P= 

14320.230 

PL= 

13207.990 

T1Z = 

974.40520 

T4Z = 

990.42670 

TM= 

970.90880 

VM= 

137.12450 

MACHM= 

.21981850 

A= 

.58344870 

DQDZ- 

12404 . 000 

Q= 

3971.8080 



Z- 

.50000000 

P= 

14163.720 

PL= 

13309.140 

T1Z= 

973.55330 

T4Z= 

989 . 69870 

TM= 

969.51760 

VM- 

147.07040 

MACHM= 

.23600390 

A~ 

. 59423330 

DQDZ= 

12500.000 

Q= 

4220.9930 



BEGIN 

SECTION # 4 





MENISCUS PRESS. DIFFERENCE 




VAPOR - LIQUID, N/SQ 

M = 

854 . 52640 



CAPILLARY LIMIT, N/SQ 

M = 

6871.4060 



Z- 

. 50000000 

P= 

14163 . 670 

PL= 

13309.140 

T1Z= 

973.55330 

T4Z= 

965.82340 

TM= 

969.51760 

VM= 

147.07040 

MACHM= 

.23600390 

A= 

. 59423330 

DQDZ= 

-5945.0640 

Q= 

4220.9930 



Z= 

. 51000000 

P= 

14191.280 

PL= 

13362. 120 

T1Z = 

973.70450 

T4Z= 

965.97160 

TM- 

969.83670 

VM= 

144.78460 

MACHM= 

. 23227340 

A= 

.58383650 

DQDZ* 

“5947.4020 

Q= 

4161 . 5360 



Z~ 

. 54000000 

P= 

14271.620 

PL= 

13516.680 

T1Z= 

974 . 14320 

T4 Z = 

966.39680 

TM= 

970.76740 

VM= 

137.97880 

MACHM= 

. 22118110 

A— 

.55109640 

DQDZ- 

“5957.7470 

Q= 

3982.9570 



Z = 

. 57000010 

P= 

14348.350 

PL= 

13664.510 

T1Z = 

974.56020 

T4Z = 

966.80070 

TM= 

971.65890 

VM= 

131.24730 

MACHM= 

.21023160 

A= 

. 51570630 

DQDZ- 

“5967.8480 

Q= 

3804 . 0720 



Z= 

. 60000000 

P= 

14421.420 

PL= 

13805.420 

T1Z = 

974.95570 

T4Z= 

967.18370 

TM= 

972.51320 

VM= 

124.58740 

MACHM- 

.19941940 

A= 

.47716730 

DQDZ= 

“5977.4390 

Q= 

3624.8870 
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PL= 

TM* 

A* 


13939.610 

973.33220 

.43506880 


Z= .63000000 P* 14491.000 
T1Z* 975.33080 T4Z* 967.54690 
VM* 117.99240 MACHM* .18873140 
DQDZ= -5986.5470 Q= 3445.4260 

Z* .66000000 P* 14557.100 PL <= 14066.880 
T1Z* 975.68550 T4Z* 967.89040 TM* 974.11910 
VM= 111.45940 MACHM* .17816200 A- .38868870 
DQDZ* -5995.1740 Q* 3265.6950 

Z— .69000010 P= 14619.830 PL= 14187.420 
T1Z= 976.02110 T4Z* 968.21530 TM* 974.87700 
VM* 104.98300 MACHM= .16770040 A* .33740190 
DQDZ= -6003.3400 Q= 3085.7160 

Z= .72000000 P= 14679.150 PL* 14301.050 
T1Z* 976.33720 T4Z= 968.52150 TM* 975.60780 
VM* 98.560410 MACHM* .15734140 A* .28011890 
DQDZ* -6011.0470 Q= 2905.4960 

.75000000 P= 14735.200 PL* 14407.960 


T1Z* 

976.63500 

T4Z= 

968.80980 

TM= 

976.31530 

VM* 

92.186580 

MACHM* 

.14707530 

A* 

.21584330 

DQDZ* 

-6018.3130 

Q= 

2725.0550 



Z= 

.78000000 

P* 

14787.960 

PL- 

14507.970 

T1Z* 

976.91440 

T4Z= 

969.08030 

TM* 

976.93710 

VM* 

85.848580 

MACHM* 

. 13689080 

A* 

.14261760 

DQDZ* 

-6025.1390 

Q= 

2544 . 3990 



Z= 

.81000010 

P~ 

14837.720 

PL* 

14601.270 

T1Z* 

977.17710 

T4Z* 

969.33470 

TM= 

977.17710 

VM* 

79.497240 

MACHM* 

.12675310 

A* 

. 5739 1190E-01 

DQDZ* 

-6031.5630 

Q= 

2363.5470 



Z- 

.84000000 

p* 

14884.210 

PL* 

14687.690 

T1Z* 

977.42190 

T4Z= 

969.57180 

TM* 

977.42190 

VM* 

73.198290 

MACHM* 

. 11669880 

A* 

-.41471180E-01 


DQDZ* -6037.5610 Q« 2182.5070 


Z* 

.87000000 

P* 

14927.530 

PL* 

14767.410 

T1Z* 

977.64960 

T4Z- 

969.79220 

TM* 

977 . 64960 

VM* 

66.942950 

MACHM* 

. 10671660 

A* 

-.15713710 

DQDZ* 

-6043.1420 

Q- 

2001.2950 



z* 

.90000000 

P* 

14967.640 

PL* 

14840.250 

T1Z= 

977.85970 

T4Z* 

969.99550 

TM* 

977.85970 

VM* 

60.728450 

MACHM* 

. 96802 010E -01 

A* 

-.29517060 

DQDZ* 

-6048.3020 

Q= 

1819.9200 



Z= 

.93000010 

P= 

15004.640 

PL* 

14906.380 

T1Z* 

978.05330 

T4Z* 

970.18290 

TM* 

978.05330 

VM* 

54.550160 

MACHM* 

. 86947200E-01 

A* 

-.46195930 

DQDZ* 

-6053.0600 

Q= 

1638 . 3990 



Z* 

.96000000 

P~ 

15038.470 

PL* 

14965.620 

T1Z* 

978 . 22980 

T4Z* 

970.35380 

TM* 

978.22980 

VM* 

48.405420 

MACHM* 

. 77147890E-01 

A* 

-.66906720 

DQDZ* 

-6057.4070 

Q= 

1456.7390 



Z= 

.99000000 

P* 

15069.240 

PL* 

15018.150 

T1Z* 

978.39010 

T4Z* 

970.50890 

TM* 

978.39010 

VM= 

42.289910 

MACHM* 

. 67 396930E— 01 

A* 

-.93133480 

DQDZ* 

-6061.3610 

Q- 

1274.9560 



Z= 

1.0200000 

P* 

15096.870 

PL* 

15063.790 

T1Z= 

978.53380 

T4Z= 

970.64800 

TM* 

978.53380 

VM* 

36.201170 

MACHM* 

. 57690200E-01 

A* 

-1.2773560 

DQDZ* 

-6064.9130 

Q- 

1093.0600 



Z* 

1.0500000 

P* 

15121.470 

PL* 

15102.710 

T1Z* 

978.66160 

T4Z* 

970.77170 

TM* 

978.66160 

VM* 

30.135310 

MACHM* 

. 48021260E-01 

A* 

-1.7500160 

DQDZ* 

-6068.0740 

Q- 

911.06570 




REVERSE VAPOR FLOW IN CONDENSER 
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z— 

1.0800000 

P= 

15142.920 

PL= 

T1Z = 

978.77280 

T4Z= 

970.87940 

TM= 

VM= 

24.089640 

MACHM= 

. 38385720E-01 

A= 

DQDZ= - 

-6070.8370 

Q= 

728.98130 


REVERSE 

VAPOR FLOW IN 

CONDENSER 


Z= 

1.1100000 

P= 

15161.250 

PL= 

T1Z = 

978.86780 

T4Z= 

970.97130 

TM= 

VM= 

18.060650 

MACHM= 

. 28777730E-01 

A= 

DQDZ= - 

-6073.2060 

Q- 

546.82100 


REVERSE 

VAPOR FLOW IN 

CONDENSER 


Z= 

1.1400000 

P= 

15176.320 

PL= 

T1Z= 

978.94580 

T4Z= 

971.04670 

TM= 

VM= 

12.045600 

MACHM= 

. 19192820E-01 

A= 

DQDZ= - 

-6075.1600 

Q= 

364.59520 


REVERSE 

VAPOR FLOW IN 

CONDENSER 


Z = 

1.1700000 

P= 

15187.950 

PL= 

T1Z= 

979.00600 

T4Z= 

971.10500 

TM= 

VM- 

6.0410900 

MACHM= 

. 96253 250E-02 

A= 

DQDZ= ■ 

-6076.6880 

Q= 

182.31780 


REVERSE 

VAPOR FLOW IN 

CONDENSER 


Z- 

1.2000000 

P= 

15194.890 

PL= 

T1Z= 

979.04190 

T4Z= 

971.13950 

TM= 

VM= 

. 45123920E-01 

MACHM= 

. 71894960E-04 

A= 

DQDZ= - 

-6077.7180 

Q= 

. 3 9939 DOOE-Ol 



WEEPING OF WICK HAS OCCURRED AT Z=* 1.1400000 

THERMAL CONVERGENCE FOUND. RESULT MAY NOT BE VALID. 


SUMMARY 

EVAPORATOR UPSTREAM SURFACE TEMP. 
EVAPORATOR UPSTREAM VAPOR TEMP. 
CONDENSER DOWNSTREAM SURFACE TEMP. 
CONDENSER DOWNSTREAM VAPOR TEMP. 
TOTAL HEAT INPUT TO PIPE, WATTS 


993.40250 

980.49690 

971.13950 

979.04190 

4220.9930 


TBOIL= 6.9875950 DELTCR 1314.1080 


LEWIS HEAT PIPE CODE, NASA LEWIS RESEARCH CENTER 
BEGIN LISTING OF INPUT DATA 


15134.730 

978.77280 

-2.4433320 


15160.040 

978.86780 

-3.5395980 


15178.440 

978.94580 

-5.5842630 


15190.120 

979.00600 

-10.523270 


15194.890 

979.04190 

********** 


CHECK LIMITS. 
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APPENDIX E 


SECOND EXAMPLE CASE 

A second example case that uses several options of the code differing from the previous case has 
been included. Figure 6 shows the boundary conditions for this case. An output from the computer 
code is also given. The input variables are listed with their appropriate values at the beginning of the 
output printout. This analysis is for a water heat pipe with an outside diameter of 4.13 cm and a 
screen wick with one artery. 

The data summary at the end of the output shows that the pipe is operating very near the boiling 
limit computed by Marcus’ equation (ref. 21). However, heat pipes of this general configuration but 
with much thicker screens have been run successfully in the laboratory. This illustrates the point that 
boiling theory, particularly with heat pipes, has not been sufficiently developed to provide accurate 
predictions of boiling onset. Experimental verification is required. 


LEWIS HEAT PIPE CODE, NASA LEWIS RESEARCH CENTER 
BEGIN LISTING OF INPUT DATA 

"G" FIELD STRENGTH, FRACTIONS OF ”G H .00000000 

INITIAL # OF RUNGE-KUTTA STEPS CHOSEN 100 
NUMBER OF DATA POINTS PRINTED ALONG PIPE 20 
NUMBER OF SECTIONS IN THE PIPE 3 


SCREEN WICK, PROPERTIES COMPUTED IN SUBROUTINE 


WIRE DIAMETER, m 
SPACE BETWEEN WICK WIRES, in 
CRIMPING FACTOR 
NUMBER OF SCREEN WRAPS 


. 80000000E-04 
. 80000000E-04 
1.0500000 
2 


WITH ARTERIES 

ARTERY EXTERNAL DIAMETER, m 
ARTERY WALL THICKNESSS, m 
NUMBER OF ARTERIES 


. 60000000E-02 
. 50000000E-05 
1 


WICKING , X-SECTION DIMENSIONS SAME AS SECTION 1 
SECTION NUMBER 1 

SCREEN WICK WITH ARTERIES SUBROUTINE 
AK= . 31599480E-10 PE= .11743750 PORE= 

AV= . 10974980E-02 DCRITE= . 16000000E-03 
ROM- . 8 00 0000 0E- 04 DO= . 41300000E-01 DI= 

RSS= .00000000 RN= . 10000000E-05 TWI= 


.58766630 

. 38500000E-01 
♦ 32000000E-03 


NA- 1 


TYPE OF PIPE MATERIAL CHOSEN 
6 = COPPER 

TYPE OF WICK MATERIAL CHOSEN 
6 « COPPER 

WEBER NUMBER 6.2831800 

WORKING FLUID 
WATER 

SECTION NUMBER 1 HEAT INPUT SPECIFIED IN THIS SECTION 

1=1 ZI= .00000000 1=2 ZI= .30000000 QSEC= 2000.0000 

SECTION NUMBER 2 HEAT INPUT SPECIFIED IN THIS SECTION 

1=1 ZI= .30000000 1=2 ZI= .40000000 QSEC= .00000000 

SECTION NUMBER 3 ENVIRONMENT TEMPERATURE SPECIFIED IN THIS SECTION 
ENVIRONMENT TEMPERATURE, K 300.00000 
ZI(1) .40000000 ZI(2) 1.0000000 

SA= .15000000 RADK= .00000000 HC= 275.00000 

UPSTREAM SURFACE T4E= 409.68140 
CALL SUMEP 
COMPLETED PASS 1 
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UPSTREAM SURFACE T4E= 409.56240 
CALL SUMEP 
COMPLETED PASS 2 

UPSTREAM SURFACE T4E= 409.54920 
CALL SUMEP 

BEGIN PRINTOUT OF RESULTS 

BEGIN SECTION # 1 
MENISCUS PRESS. DIFFERENCE 

VAPOR - LIQUID, N/SQ M = 3.8750000 

CAPILLARY LIMIT, N/SQ M = 1386.8360 


z= 

.00000000 

P= 

181301.10 

PL= 

181297.30 

T1Z= 

390.41030 

T4Z= 

409 . 54920 

TM= 

390.41030 

VM* 

.00000000 

MACHM- 

.00000000 

A= 

.59661110 

DQDZ = 

6666.6670 

Q= 

. 00000000 



Z= 

. 50000000E-01 

P- 

181301.20 

PL= 

181297.30 

T1Z- 

390.41030 

T4Z = 

409.54920 

TM= 

390.41030 

VM“ 

.13361750 

MACHM- 

. 27750460E-03 

A* 

.59661100 

DQDZ= 

6666. 6670 

Q“ 

333.33330 



Z- 

. 10000000 

p= 

181301.20 

PL* 

181297.40 

T1Z= 

390.41030 

T4Z* 

409.54920 

TM= 

390.41030 

VM* 

.26723490 

MACHM- 

. 55500900E-03 

A* 

.59661100 

DQDZ- 

6666.6670 

Q= 

666.66660 



Z= 

. 15000000 

P= 

181301.20 

PL* 

181297.50 

T1 Z- 

390.41030 

T4 Z* 

409.54920 

TM= 

390.41030 

VM* 

.40085230 

MACHM- 

. 83251340E-03 

A* 

.59661100 

DQDZ = 

6666.6670 

Q= 

999.99980 



Z= 

.20000000 

P* 

181301.20 

PL* 

181297 . 70 

T1Z* 

390.41030 

T4Z = 

409.54920 

TM= 

390.41030 

VM= 

.53446980 

MACHM* 

. 11100180E-02 

A* 

.59661100 

DQDZ* 

6666. 6670 

Q= 

1333.3330 



Z= 

.25000000 

P= 

181301.20 

PL= 

181298.00 

T1Z* 

390.41030 

T4Z* 

409.54920 

TM= 

390.41030 

VM= 

.66808740 

MACHM* 

. 13 87 52 3 0E-02 

A= 

. 59661100 

DQDZ = 

6666.6670 

Q= 

1666.6670 



Z* 

. 30000000 

P= 

181301.20 

PL* 

181298.30 

T1Z* 

390.41030 

T4Z = 

409.54920 

TM= 

390.41030 

VM= 

.80169740 

MACHM* 

. 16650120E-02 

A= 

.59661100 

DQDZ* 

6666 . 6670 

Q- 

1999.9820 



BEGIN SECTION # 2 





MENISCUS 

! PRESS. DIFFERENCE 




VAPOR 

- LIQUID, N/SQ 

M * 

2.8437500 



CAPILLARY LIMIT, N/SQ 

M * 

1386.8360 



Z= 

.30000000 

P= 

181301.10 

PL* 

181298 . 30 

T1Z= 

390.41030 

T4Z = 

390.41030 

TM* 

390.41030 

VM* 

. 80169740 

MACHM* 

. 16650120E-02 

A* 

.59661100 

DQDZ= 

. 00000000 

Q~ 

1999.9820 



Z = 

.35000000 

P= 

181301.20 

PL* 

181298.70 

T1Z = 

390.41030 

T4 Z— 

390.41030 

TM* 

390.41030 

VM- 

.80169740 

MACHM* 

. 1665012 0E-02 

A* 

.57455280 

DQDZ- 

.00000000 

Q- 

1999 . 9820 



Z= 

.40000000 

P= 

181301.20 

PL* 

181299.00 

T1Z= 

390.41030 

T4 Z* 

390.41030 

TM= 

390.41030 

VM= 

.80169740 

MACHM= 

. 1665012 0E-02 

A* 

.55347000 

DQDZ* 

. 00000000 

Q= 

1999 . 9820 
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BEGIN SECTION # 3 
MENISCUS PRESS. DIFFERENCE 

VAPOR - LIQUID, N/SQ M - 2.1250000 

CAPILLARY LIMIT, N/SQ M = 1386.8360 


z= 

.40000000 

P= 

181301.10 

PL— 

181299.00 

T1Z= 

390.41030 

T4Z- 

380.80740 

TM- 

390.41030 

VM= 

.80169740 

MACHM- 

. 16650120E-02 

A— 

.55347000 

DQDZ= 

-3333.3060 

Q= 

1999.9820 



z« 

.45000000 

P- 

181301.20 

PL- 

181299.40 

T1Z= 

390.41030 

T4Z- 

380.80740 

TM- 

390.41030 

VM= 

.73488940 

MACHM- 

. 15262610E-02 

A= 

.51142280 

DQDZ- 

-3333.3050 

Q“ 

1833.3170 



Z= 

.50000000 

P- 

181301.20 

PIr* 

181299.70 

T1Z- 

390.41030 

T4Z— 

380.80740 

TM- 

390.41030 

VM= 

.66808130 

MACHM- 

. 13875100E-02 

A- 

.45941780 

DQDZ- 

-3333.3050 

Q* 

1666.6520 



Z= 

.55000000 

P— 

181301.20 

PL- 

181300.00 

T1Z= 

390.41030 

T4Z- 

380.80740 

TM- 

390.41030 

VM— 

.60127340 

MACHM- 

. 12487600E-02 

A- 

.39382120 

DQDZ= 

-3333.3050 

Q=* 

1499.9870 



Z= 

.60000000 

P= 

181301.20 

PL- 

181300.20 

T1Z- 

390.41030 

T4Z- 

380.80740 

TM= 

390.41030 

VM= 

.53446530 

MACHM- 

. 11100090E-02 

A- 

.30909930 

DQDZ— 

-3333.3060 

Q- 

1333.3220 



Z= 

.65000000 

P“ 

181301.20 

PL- 

181300.40 

T1Z- 

390.41030 

T4Z- 

380.80740 

TM- 

390.41030 

VM- 

.46765730 

MACHM- 

. 97125790E-03 

A- 

.19644040 

DQDZ= 

-3333.3060 

Q= 

1166.6570 



Z- 

.70000000 

P* 

181301.20 

PL- 

181300.60 

T1Z= 

390.41030 

T4Z- 

380.80740 

TM- 

390.41030 

VM- 

.40084910 

MACHM- 

. 83250680E-03 

A= 

.41012880E 

DQDZ- 

-3333.3060 

Q 33 

999.99170 



Z- 

.75000000 

P« 

181301.20 

PL= 

181300.80 

T1Z— 

390.41030 

T4Z= 

380.80740 

TM- 

390.41030 

VM- 

.33404110 

MACHM- 

. 69375600E-03 

A- 

-.18400240 

DQDZ- 

-3333.3050 

Q- 

833.32670 



Z= 

.80000000 

P=- 

181301.20 

PL- 

181300.90 

T1Z— 

390.41030 

T4Z- 

380.80740 

TM- 

390.41030 

VM- 

.26723300 

MACHM- 

. 55500510E-03 

A= 

-.53202300 

DQDZ= 

-3333.3050 

Q* 

666.66160 



Z- 

.85000000 

P- 

181301.20 

PL- 

181301.00 

T1Z= 

390.41030 

T4Z- 

380.80740 

TM* 

390.41030 

VM- 

.20042490 

MACHM= 

. 41625420E-03 

A- 

-1.1253170 

DQDZ- 

-3333.3050 

Q=* 

499.99660 



Z« 

.90000000 

P- 

181301.20 

PL- 

181301.10 

T1Z= 

390.41030 

T4Z- 

380.80740 

TM- 

390.41030 

VM- 

.13361680 

MACHM- 

• 27750310E-03 

A— 

-2.3130040 

DQDZ- 

-3333.3050 

Q- 

333.33120 



REVERSE VAPOR FLOW 

IN CONDENSER 



Z= 

.95000000 

P- 

181301.20 

PL- 

181301.20 

T1Z= 

390.41030 

T4Z- 

380.80740 

TM- 

390.41030 

VM= 

• 66808610E- 

-01 MACHM= 

. 13875200E-03 

A- 

-5.6202380 

DQDZ- 

-3333.3050 

Q- 

166.66600 



REVERSE VAPOR FLOW 

IN CONDENSER 



Z- 

1.0000000 

P- 

181301.20 

PL- 

181301.20 

T1Z— 

390.41030 

T4Z= 

380.80740 

TM* 

390.41030 

VM- 

. 774937 6 0E-05 MACHM- 

. 16094360E-07 

A- 

********** 

DQDZ- 

-3333.3060 

Q* 

. 18920090E-01 
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SUMMARY 

EVAPORATOR UPSTREAM SURFACE TEMP. 
EVAPORATOR UPSTREAM VAPOR TEMP. 
CONDENSER DOWNSTREAM SURFACE TEMP. 
CONDENSER DOWNSTREAM VAPOR TEMP. 
TOTAL HEAT INPUT TO PIPE, WATTS 


409.54920 

390.41030 

380.80740 

390.41030 

1999.9820 


TBOIL= 18.944220 DELTCR 19.055030 


LEWIS HEAT PIPE CODE, NASA LEWIS RESEARCH CENTER 
BEGIN LISTING OF INPUT DATA 
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APPENDIX F 


BRIEF REVIEW OF HEAT PIPE LIMITS 

References 10, 20, and 21 contain information on heat pipe limits, from which the following 
material has been gleaned. 


Sonic Limit 

At the sonic limit the vapor flow at the end of the heat input section is choked. An approximate 
expression for the sonic limit derived by Busse is presented in reference 10. This expression neglects 
friction effects. There are two tests for choked flow in this code. One is a sign change of the deter- 
minant of the set of equations used to find the increment in the principal dependent variables p, T m , 
and A(z) for the next step of the Runge-Kutta solution. The other is the onset of an extreme negative 
temperature gradient (at choke the pressure gradient, and therefore the temperature gradient, is 
negatively infinite). These conditions occurring before the end of the heat input or adiabatic section 
cause the code to take corrective action. 


Capillary Limit 


The capillary limit is defined as the maximum pressure differential between vapor and liquid that 
the meniscus of the vapor-liquid interface can sustain before rupture. At the point where this occurs, 
liquid is no longer pumped to the surface. Dryout occurs and the pipe may overheat. The capillary 
limit is usually written as 


2a 

P '- — 
eff 


where 

p v vapor pressure 
p { liquid pressure 
o liquid surface tension 
r eff effective pore radius 

The effective pore radius governing in cases where there is a gradation of pore sizes is taken to 
be the largest in the pipe region in question. Some wick structures are reputed to be tolerant of 
meniscus failure in some pores so that local dryout does not become a problem. 

Entrainment Limit 

Entrainment occurs when vapor shear forces lift liquid out of the wick and propel it back toward 
the condenser. The Weber number is customarily taken as the index of whether or not this will occur. 
The Weber number is defined as the ratio of inertial vapor forces attempting to dislodge the liquid to 
the surface tension forces restraining the liquid. One definition of the Weber number (ref. 20) is 
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where 


We = 


pv\ 

a 


p vapor density 

v vapor velocity 

Z c dimension characterizing vapor liquid surface 

In wicked heat pipes the Weber number is related to the wick spacing. In this formulation, 
experimentation has shown that the critical Weber number at which droplets dislodge is generally 
between 2 and 3. A default value of 2 is employed in this code. The critical dimension Zg is taken in 
this code as the effective pore size of the wicking for determining the capillary limit. For 
thermosyphons and wickless heat pipes, with which this code is not concerned, different expressions 
must be used. When the entrainment limit is encountered in the code, a message is printed out and 
execution continues. It is then up to the user to alter the pipe geometry or the operating conditions 
demanded, to correct the problem. 


Viscous Limit 

The vapor flow in the heat pipe is driven by the pressure drop between the evaporator and the 
condenser. At low evaporator temperatures and low power inputs a pressure may be reached at some 
point in the condenser that corresponds to the vapor pressure at the ambient or environment tempera- 
ture. In a sense the viscous forces in the vapor have used up the pressure difference required to drive 
the vapor flow. As a result the full length of the condenser is not available for heat rejection. This 
condition is defined as the viscous limit. Because computation cannot continue, a message is printed, 
as are the data for the last completed step. In actual practice heat pipes can operate for extended 
periods at this condition before all the working fluid is lost to the inactive section. 

Boiling Limit 

Boiling in the wick occurs because of the temperature gradient established through the liquid in 
the wick when heat is being conducted. Inasmuch as the liquid-vapor interface is at saturation 
conditions, the liquid at the wall is superheated. Boiling would commence from nucleation sites 
consisting of pits or defects in the surface or the mesh openings themselves in the wick. If the sites 
are sufficiently small, surface tension enables the bubbles to resist the vapor pressure within them and 
they remain intact. At a sufficiently high temperature gradient and heat transfer rate, the bubble 
grows in diameter until it can no longer withstand the pressure. The bubble then detaches and rises to 
the surface. For determination of the boiling limit the code uses the equation presented by Marcus 
(ref. 21). The default size of the critical boiling dimension in the code is 10' 6 meters. However, if a 
screen wick or a sintered wick is chosen (wick options 2, 3, 4, and 5) the critical pore size is taken to 
be the effective pore size for determining the capillary limit. It must be pointed out that satisfactory 
information on the determination of the boiling limit is very sparse, so that a liberal factor of safety 
with respect to the limit must be used in all designs. If the boiling limit is encountered in operating 
the code, a message is printed out and computation continues. 
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Critical Point 


In some instances a set of boundary conditions may be presented for solution that will cause the 
pipe temperature to exceed the critical temperature for the working fluid. Inasmuch as the surface 
tension and the heat of vaporization vanish at the critical point, no solution is possible. Before the 
code begins execution, the initial estimate for T4E is compared with the critical temperature of the 
working fluid. If T4E exceeds an arbitrary 80 percent of critical temperature, a message is printed out 
and a new set of conditions is requested. 
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Section 



Index name 
and range 
IS * 1, NS 


Heat flow tttt^tttttl t llll 

Data ft 2 34 (See 1 23 4 5 l*1,NH<IS) 

point l 1 note) 2 1 2 3 4 5 6 


Kutta | I II III I I I llll Hill I I I I I I I I I I ISTEP 
points 


Data print I I II III III I I I I IP — 1 , 
points | | j IPSTEP + 1 

Note: One constant input, two interpolation points saved. 

Figure 1 . — Schematic of a heat pipe showing types of subdivisions 
in code. 



(a) Screen wick. 



(b) Screen wick 
with arteries. 



(c) Sintered metal. (d) Sintered metal (e) Longitudinal 
with arteries. grooves. 

Figure 2. — Heat pipe wick structure options. 
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Figure 3.— Flow chart for NASA Lewis heat pipe code. 



Figure 4. — Control volume for integral energy equation. 
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Figure 5. — Input heat loads for first example case. (Dimensions 
are in meters). 
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Figure 6.— Input heat loads for second example case. (Dimensions 
are in meters). 
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